Thermal bound entanglement in macroscopic systems and area laws 



OO 
O 

o 

(N 

Oh- 
< 

o 

m 



X3 

H— > 



00 

O 

OO 

O 



X 



Alessandro Ferraro 1 , Daniel Cavalcanti 1 , Artur Garci'a-Saez 1 , and Antonio Acm 1,2 
1 ICFO-Institut de Ciencies Fotoniques, Mediterranean Technology Park, 08860 Castelldefels (Barcelona), Spain 
2 ICREA-Institucid Catalana de Recerca i Estudis Avangats, Lluis Companys 23, 08010 Barcelona, Spain 

Does bound entanglement naturally appear in quantum many-body systems? We address this 
question by showing the existence of bound-entangled thermal states for harmonic oscillator systems 
consisting of an arbitrary number of particles. By explicit calculations of the negativity for different 
partitions, we find a range of temperatures for which no entanglement can be distilled by means of 
local operations, despite the system being globally entangled. We offer an interpretation of this result 
in terms of entanglement-area laws, typical of these systems. Finally, we discuss generalizations of 
this result to other systems, including spin chains. 
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Introduction.- Many quantum information tasks are 
based on several distant observers sharing entangled pure 
states. However, in practical situations, environmental 
noise is always present and then we unavoidably deal 
with partially entangled mixed states. The first attempt 
to overcome this degradation led to the idea of entan- 
glement distillation [3, 0| , a sequence of local operations 
assisted by classical communication (LOCC) capable of 
extracting entangled pure states from a large set of en- 
tangled mixed states. However, in some cases, the noise 
degradation happens to be irreparable: given a system 
composed by n parties there exist entangled states that 
need pure-state entanglement in order to be generated, 
nonetheless no pure-state entanglement can be recovered 
back by the n parties via LOCC. These states are known 
as bound entangled states [3] and they represent arguably 
the most striking manifestation of irreversibility in the 
context of quantum information science. 

Although several examples of bound entangled states 
have been found Q, to our knowledge they mainly have 
been inspired by mathematical intuitions and simple 
recipes to construct them are still lacking. This poses 
doubts about whether bound entangled states, despite 
their remarkable properties, are mainly a mathemati- 
cal construction. Then up to now, a question remained 
open: do bound entangled states appear inherently in na- 
ture? In particular, do these states emerge in standard 
quantum many-body systems? By standard systems we 
specifically mean local interacting systems in a macro- 
scopic thermal state characterized by only few parame- 
ters. A first attempt to address this question has been 
given in Ref . Q , where bound entanglement was detected 
in the thermal state of spin systems of up to 9 qubits. 

In this work we show that bound entangled states nat- 
urally appear as thermal states of many-body systems 
composed by a macroscopic number of particles. The 
fundamental intuition behind our results comes from the 
entanglement-area law, a property satisfied by many con- 
densed matter systems according to which the entangle- 
ment between two regions scales as the surface separat- 
ing them By relating bound entanglement to area 



laws, we provide the first construction of bound entan- 
gled states with a clear physical inspiration. In this sense, 
the route we take goes in the reverse direction with re- 
spect to the one usually pursued by many recent works 
[7|: we borrow concepts from the condensed matter field 
to get new insight on quantum information science. 

Bound entanglement and area law.- Consider a quan- 
tum system of n particles described by a local Hamilto- 
nian. For the sake of simplicity, we will restrict our anal- 
ysis to translational invariant one-dimensional systems 
of n particles. If the system obeys an entanglement-area 
law, the ground-state entanglement for a bipartite split- 
ting of the parties into two groups, say A and B, scales 
as the number of connections between them. This be- 
havior has been observed in many non-critical quantum 
systems, while logarithmic corrections may appear in the 
critical case Q. Consider now two different partitions of 
a system, one in which a contiguous half of the particles 
belongs to A and the other half to B (we will refer to such 
kind of partition as half-half), and another partition in 
which the particles with even label belong to A and the 
others to B (even-odd partition). Because of the area 
law, the entanglement will saturate for sufficiently large 
n for the half-half partition, while it will increase with 
n for the even-odd partition. In this configuration, it 
is reasonable to expect that, by increasing the tempera- 
ture, the entanglement in the even-odd partition is more 
robust to thermal noise than in the half-half partition, 
and that this behavior is preserved for large systems. 

A fundamental result in the understanding of bound 
entanglement has been to recognize that all distillable 
entangled states have a non-positive partial transposition 
(NPPT) Q . Thus if one finds a non-separable state with 
positive partial transposition (PPT) it must be bound 
entangled Now, denote by T$£ t (T|:° ) the threshold 
temperatures at which the partial transposition with re- 
spect to all half-half (even-odd) partitions becomes posi- 
tive [§] . Because of the area law, one can expect that T^ t 
is strictly smaller than T^° t . Thus, it emerges a range of 
temperatures for which the system is still entangled (as 
detected by the entanglement in the even-odd partition). 
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nevertheless single particles cannot distill pure entangle- 
ment (as the half-half partitions become PPT). Indeed, 
for any pair of particles, there is always a half-half par- 
tition for which they are in opposite sides and the par- 
tial transposition is positive according to this splitting. 
In other words, bound entangled states should appear 
in general under these conditions. The rest of the pa- 
per is devoted to put on solid grounds this intuition. We 
first consider systems of coupled harmonic oscillators and 
identify a temperature range for which bound entangle- 
ment is present. This behavior is then proven to persist 
in the macroscopic limit. We also performed the same 
calculations for spin-i models and found similar results. 

Harmonic oscillators.- Consider a system composed of 
n harmonic oscillators, each one associated with position 
and momentum operators Xi and pi respectively (i — 
1, . . . , n), described by the Hamiltonian 



XiVi. 



(1) 



The diagonal elements of matrix V describe the potential 
energy in each oscillator, while the non-diagonal terms 
give the coupling between oscillators i and j. In this sce- 
nario both the ground and the thermal states are Gaus- 
sian. In what follows the entanglement will be measured 
by the log- negativity En, which quantifies by how much 
the partial transpose with respect to a given partition 
fails to be positive [1(3]. Thus, when En = the consid- 
ered partition is PPT. In Ref. ll| En between two com- 
plementary groups of oscillators, A and B, of the thermal 
state g = exp[-#/T]/Tr{exp[-iJ/T]} at temperature T 
was obtained: 

E N = ^l^{wax[l,\ k {Q)]}, (2) 

where Q = Pu~ Puo+ , w ± = W(T)~ 1 V ± ^ , and W(T) = 
l n + 2[exp(l/ 1 /2/ T ) _ ln ]-i. We denote by {A fc (Q)}J£j 
the spectrum of the matrix Q whereas P is an n x n 
diagonal matrix with the i-th entry given by 1 or — 1 de- 
pending on which group, A or B, oscillator i belongs to. 
An exact area law for the ground-state entanglement of 
this system was proven in Ref. ■ As far as for thermal 
states, the entanglement for a given bipartition is upper 
bounded by the number of connecting points [131 ]. Here, 
we mainly consider harmonic chain systems with nearest- 
neighbors interactions and periodic boundary conditions. 
The corresponding Hamiltonian ([1]) is given by a circu- 
lant potential matrix V = circ(l, — c, 0, . . . , 0, — c), with 
< c < 1/2. The system is equivalent to a chain of har- 
monic oscillators coupled with a spring-like interaction 
and is critical when c — > 1/2. 

We used Eq. ^ to compute the log-negativity for the 
even-odd and the half-half partition [9] for different tem- 
peratures and number of particles. Our calculations show 
that the entanglement follows a strict area law for non- 
zero temperatures: it increases linearly with n for the 



FIG. 1: TSl°t (solid line) and T$£ t (dashed line) as a func- 
tion of the coupling constant c for the Harmonic chain with 
nearest-neighbor interactions composed by 800 oscillators. 
Notice that the range of temperatures T^ t < T < T^° t for 
which bound entanglement is guaranteed increases when the 
system approaches the critical point c = 0.5. The intermedi- 
ate dotted line represents the upper bound to T^-J^ as given by 
Ineq. (|13[1 (e.g., we considered m = 10 and s = 3) and valid 
in the macroscopic limit: above this line the log-negativity 
in the half-half partition is zero. The threshold T^° t in the 
macroscopic limit [given by Eq. @] coincides with the one cal- 
culated numerically for 800 oscillators (see text for details). 
As a consequence, in the shaded region we can guarantee the 
presence of bound entanglement in the macroscopic limit. In- 
set: T%° t (solid line) and T^ t (dashed line) as a function 
of the number n of oscillators composing the system (log-lin 
scale), for c = 0.3 (the same behavior is found for other values 
of c). The gap T%l° t — Tjist is seen to remain constant with 
the size of the system, apart from an initial transient. 



even-odd case, while it saturates for the half-half parti- 
tion. The temperature just defines the rate the entangle- 
ment increases with n for the even-odd partition and the 
entanglement saturation value for large n for the half- 
half partition. As shown in Fig. [T] T^° st is strictly larger 
than T££ t , indicating, as discussed above, the presence 
of bound entanglement. We performed our calculations 
for systems composed by up to 800 oscillators and found 
that all the computed threshold temperatures, and so 
also the gap T^° t — T^ t , are independent of the size of 
the system, as can be seen in the inset of FigfT] 

Macroscopic systems.- The results displayed in Figfl] 
strongly suggest that bound entanglement persists in the 
macroscopic limit. Actually, we can prove this statement 
by establishing an analytical formula for En in the even- 
odd partition and an upper bound for it in the half-half 
partition, when n — * oo. We present here the main steps 
of this proof. A more detailed analysis is given in 14 1 . 

We start with the calculation of the even-odd log- 
negativity, following the one given in Ref. [ll[ for the 
ground state. Consider the Hamiltonian fl} with V de- 
scribing nearest-neighbor interactions as shown before. 
Matrix V is circulant and the matrices uj^ can be di- 
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agonalized by a discrete Fourier transformation, imple- 
mented by a matrix fl. Namely, we have that fluj^Q^ = 
D ± , with D± = S k i 4 and 4 = A* 1/2 tanh(VXT/2T), 
where A& = 1 — 2c cos (27rfc/n) are the eigenvalues of V. 
Concerning matrix P, one has that: 



ftPO T = P = 



1 



1 

n/2 



n/2 





(3) 



The spectrum of Q coincides then with the one of 
PD~ PD + , which in turn is straightforward to calculate 
due to its block diagonal structure. Defining the function 



f(k,n,c,T) = 



I A 



fe+n/2 
~A 



tanh • 



/AT 

■ tanh ■ 



z/2 



A- 



2T 



2T 



one has that the eigenvalues of Q that can contribute 
to the log-negativity © are given by f(k,n,c,T), with 
double multiplicity and fc = 0, . . . , n/A (for n multiple of 
4, n > 4). The log-negativity of the even-odd partition is 
different from zero when the temperature T is such that 
/(0, n, c, T) > 1. In particular, the curve /(0, n, c, T) = 1 
gives the threshold temperature T£l° t , in formula: 



1 + Nanh(^W^5^]=l. (4) 



1 - 2c 



2T 



2T 



The threshold above coincides with the one depicted in 
Fig. [T] (solid line) and it is independent on the total num- 
ber of particles n (see also the solid line in the inset of 
Fig. Q}, i.e. it holds in the macroscopic limit. For tem- 
peratures below T^f st there exists a k(n,c,T) such that 
f(k,n,c,T) > 1 for k < k(n,c,T), which in turn gives 
rise to the following expression for the log-negativity: 



fc(n,c,T) 



E N =22 lo §2 f( k , n , c , T ) - / da;lo g2/( a; ! c J T )) 

where, for large n, we have replaced the sum over k by an 
integral over x — 2nk/n. As a consequence we see that 
the log-negativity grows linearly with the system size also 
for non-zero temperatures. 

Regarding the half-half partition it is possible to find 
an exact upper bound for the threshold temperature 
Tdist ■ Remarkably, this allows to identify a range of tem- 
peratures for which the presence of bound entanglement 
can be guaranteed also in the macroscopic limit. We pro- 
ceed as follows. Let us define the matrix 



i/2-l n-1 



X 



ij - w ij 



(SikSjh + SjkSih) (5) 



k=0 h=n/2 



where 8ij denotes the Kronecker delta. Following 
Ref. [l3| , the log-negativity is zero when 

^min 

[W{T)}- 2 + 2max l \X i [Xuj + }\ < 1 , (6) 



where 



Xnun[W(T)} 



3 \/l+27/T , 1 



; \/T+2H/T _ I ' 



(7) 



Recognizing that the second term in the left hand side 
of (J6j) is twice the spectral radius r(Xu> + ) of the ma- 
trix Xlu + , we can use any matrix norm to bound it 
from above [HI]. An upper bound for r{Xuj + ) is then 
given by r(Xuj+) < \\Xlo+\\ < \\X\\ ||w+||. We consider 
the maximum row sum matrix as a representative norm: 
\\A\\ = maxi%2AAij\ . The goal, now, is to bound \\uj + \\ 
and ||X||. Before proceeding, recall that w ± are circu- 
lant matrices, hence completely specified by their first 



row uj = circ(t>Q 



One can show that 



i r 27r 

1 ^ Jo 



(8) 



when n — > oo. As a consequence, for any integer s, by 
integrating by parts s times we have: 



1 



1^1 ^ 2^ 



Let us bound first 



2 71 



dx 



-^d ± (x) 
dx s v ; 
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2irl s 



(9) 



||. Being lo + a circulant matrix, 
it follows that there is no need to look for the maximum 
over the rows. Then one can write, for any integer m, 
||w + || = 5+ + £+, where we defined the partial sum and 
the residual term, respectively, as follows: 



S 



+ — 



Eh 

l=—m 



E i 

Z=m+1 



(10) 



In order to obtain a bound on ||w + || one can now fix m, 
calculate explicitly 5+ and bound £ + from above. This 
latter step can be achieved by using the bound in Eq.©, 
leading to £+ < CfC{s,m + (where ((s,m + 1) is 
the generalized Riemann zeta function). Summarizing, 
for any integer s and m, we have proven that 



C 



<S+ + ^C(s,m + l) = K+ 



(11) 



the bound above being tighter for large m and s. 

Concerning the upper bound for ||A|| we have 
that ||A|| coincides with ||_B|| where B is the sym- 
metric Toeplitz matrix with the first row given by 
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,/2' V n/2-l' ■ ■ ■ ' V 2 > V l 



We then need to determine 



which row of B has the maximum sum. Fortunately, 
when n — > oo the periodic boundary conditions can be 
disregarded and the first row of B determines its norm. 
This is because going from the first to the second row 
we simply remove the term \v^\, and so on for the other 
rows. Similarly as we did for ||cl> + || we obtain that for 
any integer s and m 

\\X\\<S- + ^<;(s,m + l) = K- ia , (12) 
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FIG. 2: Negativity in the even-odd (full symbols) and half- 
half (empty symbols) partitions for thermal states of n spin- 
one-half particles with nearest-neighbor Hamiltonian given by 
Hxx = - E(«+i+"Wti) +B£<7f, and B = 1.9. The 
temperatures for each partition are set to be, from top to bot- 
tom, T — 2, 2.6 . In the even-odd partition we can clearly see 
an increase of the negativity with respect to the system size, 
whereas it saturates in the half-half partition. Both behaviors 
are expected from the entanglement-area law. For T = 2.6 we 
see that the state is bound entangled (the negativity in the 
half-half partition being zero). 



where S~ = Y^iLi \ V T\- ^ s done before, S~ can be cal- 
culated explicitly and the bound is tighter for large m 
and s. Summarizing, considering Eqs. ©, (J7J), (fTTj) . and 
(|12[) . we have shown that, in the macroscopic limit, the 
log-negativity in the half-half partition is zero when the 
following inequality is satisfied: 



1K + K~ 

m,s m,s 



< 1 



(13) 



Based on the formula above and on the threshold Q , we 
depicted in Fig. [1] the region in the c — T plane for which 
bound entanglement is present in the macroscopic limit 
(shaded region). We see that for any coupling we can 
guarantee that there is a range of temperatures for which 
the log-negativity in the half-half partition is zero, nev- 
ertheless the state is entangled. Thus the temperature, a 
single macroscopic and measurable quantity, clearly de- 
termines the distillability properties of the system. 

Concluding remarks.- We remark that our results, 
when combined with those of Ref. [lj|, prove that en- 
tangled states such that all bipartitions are PPT, cannot 
be obtained for the harmonic systems studied here [l7j |. 
This is because the threshold temperatures for the even- 
odd partition given above coincide with those found in 
Ref. [16[ for full separability. Thus, thermal states be- 
come PPT and fully separable at the same temperature. 

Finally, we have considered other models of harmonic 
chains as well as spin systems. All the obtained results 
are consistent with the previous reasoning 14}: there is 
a temperature range for which the negativity in the half- 



half partitions is zero, nevertheless the system is still en- 
tangled as proven by the even-odd negativity. For spin 
systems we could not go beyond 12 particles due to com- 
putational hardness, but the observed numerical results 
again support the existence of thermal-state bound en- 
tanglement in the macroscopic limit (see Fig. [5]). Due to 
the generality of the area laws, the present results are ex- 
pected to hold for a variety of other systems. These find- 
ings, besides showing novel aspects of entanglement in 
many-body systems, are of relevance from an application- 
oriented perspective. As shown here, quantum correla- 
tions emerge "gradually" in a many-body system: start- 
ing at high temperatures in a classically correlated sce- 
nario, the system passes through a regime where quan- 
tum correlations start to be present, but in a bound form. 
Then, for lower temperatures, the system enters a truly 
quantum correlated regime where free distillable entan- 
glement is present and potentially available for quantum 
information applications. 

Note added.- By completing this work, bound entan- 
glement was independently found in three-qubit reduced 
states of the XY model at nonzero temperature [HI]. 
Here we have shown the existence of bound entanglement 
in the whole thermal state of the macroscopic system. 
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